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ABSTRACT 

We use the SWIFT Burst Alert Telescope (BAT) sample of hard x-ray selected 
active galactic nuclei (AGN) with a median redshift of 0.03 and the 2MASS J 
and K band photometry to examine the correlation of hard x-ray emission to 
Eddington ratio as well as the relationship of the J and K band nuclear luminosity 
to the hard x-ray luminosity. The BAT sample is almost unbiased by the effects 
of obscuration and thus offers the first large unbiased sample for the examination 
of correlations between different wavelength bands. We find that the near-IR 
nuclear J and K band luminosity is related to the BAT (14 - 195 keV) luminosity 
over a factor of 10 3 in luminosity ( L IR « L} 3 2 f T )and thus is unlikely to be clue 
to dust. We also find that the Eddington ratio is proportional to the x-ray 
luminosity. This new result should be a strong constraint on models of the 
formation of the broad band continuum. 

Subject headings: galaxies: Seyfert, surveys, X-rays: galaxies 

1. Introduction 

The estimate of black hole masses and the ratio of luminosity to the luminosity de- 
rived from the black hole mass (the Eddington ratio) are the fundamental parameters in 
the evolution of active galaxies. The Eddington ratio is defined as L/L Edd , where L Edd = 
( Mbh/Mq ) x 1.3 x 10 38 er<?s _1 or the limit where gravitational and radiation pressure are 
balanced. It is important because it is a measurement of the ratio of energy output of a 
source to the physical limit of energy output. Understanding the energy output of active 
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galactic nuclei (AGN) is essential to understanding both the nature of AGN sources and the 
effect they have on their environment. Therefore, measuring both and the Eddington 
ratio are the subject of much recent work (Kollmcier et al. 2006; Netzer & Trakhtenbrot 
2007; Woo & Urry 2002), based primarily on optical samples and the combination of x-ray 
and IR samples. These papers have either used data from the literature (Woo & Urry 2002), 
the SDSS (Netzer & Trakhtenbrot 2007), or deeper optical, x-ray and IR data (Kollmcier et 
al. 2006) to obtain conclusions about the relationship of the luminosity to black hole mass 
and its evolution across cosmic time. For example, Kollmeier et al. (2006) conclude that 
there is strong evidence that supermassive BHs gain most of their mass while radiating close 
to the Eddington limit, and that the fueling rates in luminous AGNs are determined by BH 
self-regulation of the accretion flow rather than galactic-scale dynamical disturbances. 

Whereas previous samples are strongly biased against the most numerous AGN in the 
universe, those with high line of sight column densities, e.g. ‘obscured’ AGN (Matt et al. 
2006), our BAT sample, since it is selected solely based on 14-195 keV flux, is not. However 
the BAT sample, as opposed to most optical samples, has a very low median redshift (0.03) 
and very few objects at z > 0.1. Since there is a strong relationship between line of sight 
obscuration and intrinsic luminosity (Tozzi et al. 2006; Markwardt et al. 2005), such that the 
less luminous sources are much more likely to be obscured, the optical samples, which miss 
most of the obscured sources, are strongly deficient in low luminosity sources. Furthermore, 
there is often a complex relationship between measures of obscuration in the x-ray, IR, and 
optical bands and there may be strong evolutionary effects such that the fraction of obscured 
objects is higher at higher redshifts (La Franca et al. 2005). If there is a physical relationship 
between luminosity or stage in the life of the AGN and the occurrence of obscuration, as 
predicted by several models (Hopkins et al. 2006), then the results from such surveys may be 
rather different than The BAT sample which is heavily weighted to the local universe with 
a mean redshift of « 0.03. 

In addition, there has been no conclusive determination of the origin of the ‘non-thermal’ 
continuum in AGN. In the x-ray band, there are two mechanisms that have received exten- 
sive discussion: thermal Comptonization of photons produced by the accretion disk from 
a population of high energy particles, perhaps in a corona above the disk (Zdziarski et al. 
1990), and emission due to jet processes (Kording, Falcke, & Corbel 2006). However, there 
is at present no agreement as to the origin of ‘non-thermal’ IR and optical emission in radio 
quiet objects or even whether it exists. Most models for the origin of the optical emission 
have focused on disk models (e.g. Sun & Malkan (1989); Czerny & Elvis (1987)) with some 
Comptonization. Dusty torus models (Pier & Krolik 1993), as the origin of the K band 
flux, have received support from reverberation observations (Suganuma et al. 2006) in the 
optical (UBV) and near-infrared (JHK) bands. These studies found a clear time-delayed 
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behavior of the K band with respect to the V band, strongly suggesting that the bulk of 
the K-band emission originates from thermal emission of dust grains in an optically thick 
torus (Glass 2004). There are strong correlations between the optical and x-ray properties 
(Boroson 2002; Brocksopp et al. 2006) and the IR properties (Alonso-Herrero et al. 2001), 
whose origin is not easy to explain if these continuum components originate from completely 
different physical mechanisms as is predicted by the most modern dusty torus models of the 
IR continuum (Nenkova et al. 2008). 

In this paper we present the first results on the Eddington ratio from an unbiased survey 
of AGN selected in the 14 - 195 keV band with the Swift BAT survey. We find a strong 
correlation of the Eddington ratio to the hard x-ray luminosity over a wide range in L/L em 
(10“ 4 to 3) in our sample, as opposed to the narrow range in the sample of Kollmeier et 
al. (2006) selected from the AGES survey. In addition we find that the luminosity of the 
point-like component in the 2MASS J and K band data is almost linearly correlated with 
the high energy luminosity indicating that these two components are very closely coupled 
and probably have the same physical origin. Since the x-ray emission is not due to dust 
processes, these results indicate most of the near IR is not either. 


2. Survey 

The BAT 9 month survey (Tueller et al. 2008) is an all sky survey in the 14 - 195 keV 
band (Markwardt et al. 2005) with a sensitivity threshold that varies from 1— 3x 10 -11 erg cm -2 s^ 1 
in the 14 - 195 keV band, as a function of sky exposure. At \b\ > 15°, the survey is almost 
complete with 130 total sources of which 103 are identified AGN and the rest are galac- 
tic objects with only 1 unidentified source. Of the 103 AGN, 14 are Blazars and will not 
be included in the following discussion. All but 3 of the identified objects have redshifts. 

The survey is continuing and the sensitivity is increasing roughly as the square root of time 
(Tueller et al. 2008). 


3. Black Hole Mass Estimation 

While there are several published papers on scaling black hole masses to optical and 
IR luminosities, we have chosen to use the recent paper by Novak, Faber, & Dekcl (2006). 
They give scaling relations between the central black hole mass and global galactic properties 
based on optical and IR data. Since, for our sample, the only uniform optical or IR data 
set is the 2MASS data and we wish to minimize the contribution of possible AGN light to 
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the total IR light, we use J and K band magnitudes. As is well known, the ratio of the 
host galaxy to the QSO photon flux is at a maximum in the NIR, because the SED of many 
stellar populations has a maximum at 1 . 6 /jrn and the AGN flux has a minimum at 1.2/mi 
(Elvis et al. 1994; Richards et al. 2006). The AGN flux increases at shorter wavelengths due 
to accretion-disk or power-law emission and at longer wavelengths due to thermal emission 
of hot dust heated by the AGN and/or starformation contributions.. 

We obtained the total J and K 2MASS magnitudes (m j(tot) and m K {tot)) for the BAT 
AGN sample, as listed in NED. Additionally, we used the IRAF task qphot to obtain circular 
aperture photometry of the central region, equal to the PSF FWHM of the 2MASS images. 
While we label these values as the nuclear magnitudes (m j(nuc) and m k(uuc)), they are 
clearly overestimates to the nuclear flux, since no provision is made for the contribution 
of the galaxy light in the central region. We have found that in most objects the possible 
contribution to the total J band light from the nucleus is small and in most cases is well 
measured. This is important, since if the J band nuclear light dominated the emission, our 
estimates of BH mass would be severely impacted and the scatter in relations involving the 
BH mass (such as the Eddington ratio) would have much more noise. 

From the apparent magnitudes, we computed absolute magnitudes (M j(tot), M xitot), 
M j(nuc), and M x(nuc)) with the redshifts listed in Tuellcr et al. (2008) and assuming H 0 = 
75kms -1 Mpc -1 . From the total and nuclear absolute magnitudes, we computed estimates 
of the galactic/stellar component of the NIR emission (M j (stellar) and M x (stellar)) as 
Msteiiar = 2.51og(x/(x - 1)) + M tot where x = F tot /F nuc = 10 -o.4(M tot -M nuc ). Following 
Novak, Faber, & Dekel (2006), we obtain black hole masses for each AGN in our sample 
using their relation: log (M(BH)) = 8.19 + 0.524 x (—M K ^ ste u ar ) —23) or that logM(£>77) « 
1.3 log L(K ste u ar ). Our estimated black hole masses can be found in Winter et al. (2008), 
along with the 0.3-10 keV X-ray properties. 

We note that aperture photometry in the J and K band (Ward et al. 1988; Glass & 
Moorwood 1985) allows one to directly compare the AGN J,K band fluxes to the stellar flux 
for a sub-sample of our objects. We find that, in general, the total (stellar +nuclear) J and K 
band flux is dominated by star light with the AGN contribution, even for the most luminous 
BAT selected AGN, being less than « 30% of the total luminosity. This dominance of the 
near IR light by stars is beautifully confirmed by the spectroscopic IR data for PG quasars 
(Dasyra et al. 2007). Even in these very luminous quasars (all of which have non- stellar 
IR luminosities more than 3 times larger than the stellar luminosities of the most luminous 
object in our sample and 80 times larger than most of the objects), about 10% of the total H 
band light is due to stars. The work of Bentz et al. (2006), which contains some of the same 
objects as in our sample, shows that even in the optical band at 5100A (where the relative 
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dominance of AGN light is larger than in the near IR) none of their objects are dominated 
by AGN light. We thus believe that the noise introduced by the measurement error of the 
AGN light on the estimate of black hole masses for these objects is relatively small (Peng et 
al. 2006). We have compared our photometry with that in the Peng et al. paper for both the 
extended emission and the nuclear emission. Peng et al. perform a detailed deconvolution 
of the 2MASS images and thus, presumably have done a more careful job than we have. 
However we find only « 1/3 mag scatter in the K band nuclear magnitudes between the 
Peng results and ours, indicating that our technique is sufficient to the task at hand and our 
relatively coarse photometry does not effect any of the results in this paper. 


4. Bolometric Correction 

We use a constant ratio of 15 to transfer the Swift BAT (14 - 195 keV) luminosity 
to bolometric luminosity, based on the broad band SEDS (Barger et al. 2005) and the 
measurement of an E -2 photon index for the Swift BAT sources over the entire 14 - 195 keV 
band. While there seems to be a relation between the softer band x-ray luminosity and 
bolometric luminosity (Steffen et al. 2006), the effect is rather shallow and we choose to 
neglect it in this paper since it has not yet been measured in our harder x-ray band. If the 
effect is similar in the BAT band to that in the 0.5-2 keV band, it will tend to increase our 
estimate of the bolometric luminosity of the higher luminosity objects. We are accumulating 
broad band SEDS for all of these objects, but for many of the objects the nucleus is very weak 
in the UV and optical bands and thus measuring its luminosity is difficult. To obtain good 
quality estimates of the nuclear optical and UV fluxes will probably require HST (Malkan 
et al. 1998) or adaptive optics data. 


5. Results 
5.1. Eddington Ratio 

Using the apparent magnitudes defined in Section 3, we derived associated NIR fluxes 
with the flux for zero- magnitude zero point conversion values of 3.129 x 10 -10 ergs -1 cm -2 
for the J band and 4.283 x 10~ n ergs -1 cm -2 for the Kg band (Cohen et al. 2003). Luminosi- 
ties, for example the nuclear luminosity L j( nuc ) and the galactic/stellar luminosity L K (stellar), 
were computed, as for the absolute magnitudes. Additionally, we computed the bolometric 
luminosity (bolometric correction xLbat) and the Eddington luminosity (using the esti- 
mated black hole mass from Mr (stellar) and the Novak, Faber, & Dekel (2006) equation) to 
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Fig. 1. — Plots of (a) the J filter nuclear luminosity versus the BAT x-ray luminosity (14 - 
195 keV), (b) the K filter stellar luminosity versus the BAT x-ray luminosity, (c) the ratio 
of the bolometric luminosity (estimated from the BAT luminosity) to Eddington luminosity 
(estimated from M(BH) (Novak, Faber, & Dekel 2006) using the K filter stellar (total - 
nuclear) magnitude) (L/L EM ) versus the mass derived from the stellar K filter magnitude 
(in M 0 ), and (d) the ratio L/L em versus the BAT luminosity. The circles represent Syl - 
Syl.2 sources, the squares represent Syl.5-Syl.9, and the triangles represent Sy2s. For most 
of the plotted points, the error bars fit entirely within the size of the symbol. The lines in 
panels (a) and (d) represent linear fits to the data using the OLS bisector method (results 
are discussed in the text). 
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determine the Eddington ratio for each source. 

We show in Figure 1 (a) the correlation of BAT x-ray luminosity to L j( nuc ), (b) the 
relation of BAT x-ray luminosity to L k ( stellar), ( c ) the Eddington ratio versus black hole mass 
(M(BH)/Mq), and (d) the correlation of the Eddington ratio to the BAT x-ray luminosity. Of 
the final 3 panels (b,c,d in Figure 1), only the Eddington ratio to the x-ray luminosity shows a 
significant correlation, of the form logL bo ; /L Edd — (0.924±0.069) logL Bj 4 T +(— 42.055±2.999) 
(solid line in Id) with a correlation coefficient of R2 = 0.45 for 98 points which has a 
probability of 3x 10~ 7 by chance. Notice that both absorbed and unabsorbed AGN (primarily 
type II and I respectively) fall on the same scatter diagram indicating the absence of Compton 
thick objects, which might be expected to have a high IR to x-ray ratio. Our results are 
similar to previous work (Woo & Urry 2002) where no correlation of Eddington ratio to 
black hole mass is found. However, our median Eddington ratio of 0.069 is much less than 
that found by Kollmcier et al. (2006) who find values of ~ 0.2. We attribute this to the 
much higher luminosity of the Kollmcier et al. (2006) sources, since we find for the most 
luminous objects in our sample a similar Eddington ratio of 0.2. The primary relationship is 
clearly with the BAT hard x-ray luminosity, since our assumption of a constant conversion 
from BAT to bolometric luminosity can only add noise, as does our method of estimating 
the black hole mass. 

The prime systematic errors are our estimate of the black hole mass and the bolometric 
correction. These should, on average, not reduce the significance of the relationship between 
the Eddington ratio and the hard x-ray luminosity. The sense of the errors is to have 
overestimated the black hole mass clue to the improper correction for AGN light (e.g. the 
black hole mass is maximized in using the Nowak correlations if more of the K or J band 
light is assumed to be clue to stars). This error will be a bigger problem for the higher 
luminosity objects because more of the K band light in the 2MASS aperture is clue to the 
AGN. Since the hard x-ray luminosity does not have a similar possible error, we may have 
systematically underestimated the black hole mass and thus overestimated the Eddington 
ratio in the high luminosity objects. Similarly, if in the BAT band the x-ray light is a smaller 
fraction of the bolometric luminosity for luminous objects, as it seems to be in the 0.5-2 keV 
band for type I objects (Steffen et al. 2006), then the true bolometric luminosity is larger at 
the higher L(BAT). Both of these errors will tend to reinforce the sense of the correlation. 
Alternatively, if we have overestimated the AGN J band light it will reduce the total IR 
stellar light and thus our mass estimates and effect (see next section) the relationship of IR 
nuclear to x-ray light. Since this relation shows a scatter of only « 0.5 dex, this effect does 
not appear important. Of course there is some concern (as pointed out by the referee) that 
this correlation is an artifact of relating the x-ray luminosity to the x-ray luminosity divided 
by the black hole mass, where the black hole mass has a small range and the x-ray has a 



large range. We do not believe that this is the origin of our correlation since the correlation 
of L(BAT) with stellar mass has the same narrow range in mass and large range in L(x) 
and there is no correlation. Almost all correlations with black hole mass have the same bias, 
since there is only an approximately factor of 300 range in black hole masses estimated from 
bulge luminosities, one should exercise caution. 


5.2. Correlation of J and K band nuclei luminosities with hard x-ray 

As is obvious from Figure 1, the luminosity of the nucleus in the BAT and J band 
is strongly correlated, with a best fit (using the ordinary least square bisector method as 
presented in Isobe, Feigelson, & Nelson (1986)) of log (Lj (ram)) = log (Lbat) x (1.257 ± 
0.077) + (—6.57 ± 3.09) and R 2 = 0.63. The same correlation is seen with the K band 
nuclear luminosity (not shown). Using the same OLS bisector method, Steffen et al. (2006) 
finds a correlation between optical (2500A) and soft x-ray (2 keV) luminosity with a slope of 
0.721 ±0.011 in a sample of unobscured, optically-selected AGN. The non-linear relationship 
between optical and soft x-ray emission is partially true due to the fact that low x-ray 
luminosity sources are more likely to be absorbed (Steffen et al. 2003). Since the very hard 
x-ray /BAT and IR luminosities are relatively insensitive to obscuration, the almost linear 
relationship between these luminosities is not sensitive to obscuration effects and is thus 
more robust than the optical or UV to soft x-ray correlations. 


6. Discussion 

The median value of the Eddington ratio that we find is L/L Edd ~ 0.02 with a full range 
of « 1000. This value seems to be at odds with the work of recent AGN samples, such as 
the Kollmeier et al. (2006) sample selected based on optical/IR/x-ray properties, which finds 
Lboi/ LecLcL = 0.25 over a rather different luminosity range 10 45 - 10 1 ' ergs -1 in L(Bol) than 
the BAT sample. We believe that this is because the Kollmeier selection is really an optical 
one, since only the broad line sources have mass estimates. Such a sample is heavily biased 
towards “unobscured” sources which have higher luminosities (Steffen et al. 2003). Also, we 
now know that they also have higher L/L^dd (Winter et al. 2008). Additionally, the Kollmeier 
selection is also biased towards much higher redshifts than the BAT survey, which could also 
explain the difference if, for instance, their is an evolution in Eddington ratios. The wide 
range of inferred Eddington ratios from x-ray data has been noticed before (Merloni, Heinz, 
& di Matteo 2005) and is also seen in other low z samples (Hopkins, Narayan, & Hernquist 
2005). However, our sample differs significantly from the Ho, Filippenko, & Sargent (1995) 
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sample in not having objects with very low Eddington ratios. We are not sure why we 
would have a bias towards higher Eddington ratios than Ho, Filippenko, & Sargent (1995). 
However, since we have not convolved our selection function with a theoretical Eddington 
ratio model, one cannot be sure that our differences with Ho and Kollmcir are not due to 
selection biases, the luminosity range covered, and/or evolutionary effects. 

Softer x-ray surveys in the 2-10 keV band (e.g. Steffen et al. (2006)) have found that 
the bolometric correction is a function of luminosity, such that the x-ray band carries less of 
the luminosity at higher luminosities. If this were also true for our sample it would increase 
the amplitude of the correlation between hard x-ray and bolometric luminosities and change 
the slope (Hopkins, Narayan, & Hernquist 2005). 

The correlation between Eddington ratio and x-ray luminosity has, to our knowledge, 
not been reported from other samples (e.g. Merloni, Heinz, & di Matteo (2005)) but it is 
clearly present in them. If one graphs the Merloni, Heinz, & di Matteo (2005) data one 
finds exactly the same correlation as in our plot. There are theoretical reasons (Hopkins, 
Narayan, & Hernquist 2005) to believe that the most luminous objects have the highest 
Eddington ratios and the lower luminosity objects have lower L em- The Hopkins paper 
gives the Eddington ratio as a conditional probability distribution in which the width is a 
function of Eddington ratio - this is hard to interpret in our present framework and we have 
not tested the predictions of this model. 

We note that our distribution of masses is weighted towards higher masses than predicted 
by Hopkins, Narayan, & Hernquist (2005) at low z, but we cannot conclude that this is 
a strong result since we have not convolved the Hopkins, Narayan, & Hernquist (2005) 
prediction with our sensitivity limit. 

It is well known that there is a break in the luminosity function of AGN, the origin of 
which is not clear. Below the break in the luminosity function, the median Eddington ratio 
from our sample is -2.1 and the variance is 0.35, while at higher L the median is -1.64 and 
the variance is 0.199. This is in agreement with the predictions of Hopkins, Narayan, & 
Hernquist (2005), that the variance in the Eddington ratio should decrease at higher L, and 
provides a possible explanation of the break in the luminosity function as a change in the 
mean Eddington ratio of the population. 

Our finding of a strong almost linear correlation of the near IR luminosities with the 
BAT luminosity (see Horst, et al. (2007) for a similar result) argues that similar processes 
produce the near IR continuum and the hard x-rays. Since the mechanism producing the 
x-ray emission is unrelated to dust, these results do not favor dust reprocessing models for 
the origin of the near IR continuum, unless the main driver of the dust reprocessing is simply 
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and linearly connected to the x-ray continuum intensity. 
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